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ABSTRACT 

The Hyades K2V+WD system 471 Tau is a prototype post-common envelope 
system and a likely cataclysmic binary progenitor. We present 10 days of nearly 
continuous optical photometry by the MOST (Microvariability & Oscillations of 
STars) satellite and partly simultaneous optical spectroscopy from DDO (David 
Dunlap Observatory) of the binary. The photometric data indicate that the spot 
coverage of the K dwarf component was less than observed in the past, suggesting 
that we monitored the star close to a minimum in its activity cycle. Despite the 
low spot activity, we still detected seven flare-like events whose estimated energies 
are among the highest ever observed in V471 Tau and whose times of occurrence 
do not correlate with the binary orbital phase. A detailed O — C analysis of the 
times of eclipse over the last ~ 35 years reveals timing variations which could be 
explained in several ways, including perturbations by an as-yet-undetected third 
body in the system or by a small orbital eccentricity inducing slow apsidal motion. 
The DDO spectra result in improved determinations of the K dwarf projected 
rotation velocity, V/^sinz = 92 km s~^, and the orbital amplitude, = 150.5 
km s~^. The spectra also allow us to measure changes in Ha emission strength 
and radial velocity (RV) variations. We measure a larger Ha velocity amplitude 
than found previously suggesting that the source of the emission in V471 Tau 
was less concentrated around the sub-white-dwarf point on the K star than had 
been observed in previous studies. 
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Subject headings: stars: close binaries - stars: eclipsing binaries - stars: variable 
stars - photometry: space based 



1. INTRODUCTION 



V471 Tau is a close eclipsing binary star {V ~ 9) consisting of a hot white dwarf 
and a red dwarf with a n orbital period of 0.521 d. It is a member of the Hyades cluster 
( Werner fc Rauchlll997l ) and - very likely - a cataclysmic binary progenitor (IStill fc Hussain 
20031 ). V471 Tau has bee n the subject of numerous in v estigations over the p a st 35 years; cf. 



the m ai n contributions bvlSkillman fc Patterson 



( I2OOII ) , llbanoglu et al.l (120051 ) and iHussain et al. 



1988h.lciemens et al.Nl992h . lO'Brien et al.l 



(120061 ). 



The V471 Tau system may be the prototype of a post-common envelope binary with 
a white dwarf and a main sequence star. The mass and radius of both components can be 
measured with high accuracy, while the K dwarf which is spun up to high rotation rates by 
tidal forces may be an analogue for rapidly rotating pre-ZAMS stars like AB Dor. Simultane- 
ous precise time-resolved photometry and spectroscopy which cover phases of eclipse in the 
V471 Tau system can sample the spot coverage of the K dwarf. Eclipse timing can measure 
apsidal motion in the binary and test whether the system is actually a triple one with a 
third undetected component. We therefore organized a coordinated campaign of spacebased 
photometry from the MOST satellite and groundbased spectroscopy from DDO. 

We present new MOST and DDO observations of V471 Tau in Section d The MOST 
light curve and its changes are discussed in Section [3] while Section H] gives a description of 
the spectroscopic observations. Section [5] summarizes the combined results. 



^Based on data from the MOST satellite, a Canadian Space Agency mission jointly operated by Dynacon 
Inc., the University of Toronto Institute for Aerospace Studies and the University of British Columbia, with 
the assistance of the University of Vienna, and on data obtained at the David Dunlap Observatory, University 
of Toronto. 
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OBSERVATIONS OF V471 TAU 



2.1. MOST photometry 



The MOST (Micr ovariability & OsciUations of STars) space mission (jWalker et al.ll2003l : 
Matthews et al.ll2004j ) was designed to perform high-precision optical photometry of bright 
stars with long time coverage and high duty cycle. MOST is equipped with a 15-cm telescope 
and a custom broadband filter (spectral transmission peak ~ 5500 A and FWHM ~ 3000 A). 
The polar Sun-synchronous orbit of the satellite allows it to monitor stars in the anti-solar 
direction for up to 60 days without interruption. 

MOST observed V471 Tau for 10.0 days during 4-14 December 2005 (in Terrestial Time 
Julian Days: 2,453,708.5117 - 2,453,718.5122, see below in Section [321) , covering just over 
19 orbital peri ods of the binary system. The photometry was obtained in MOST's Direct 
Imaging mode (IRowe et al.ll2006l ). with a slightly defocused stellar image sampled in a CCD 
sub-raster. The exposure time was 6.52 s, sampled at 10-s intervals. Two reference stars in 
the same field (CSC 01252-00692, V = 8.9 and CSC 01252-00046, V = 9.8) were observed 
simultaneously in the same way to calibrate instrumental or satellite orbital artifacts. 

The MOST instrument focal plane can be illuminated by scattered Earth light whose 
level is modulated by the MOST orbital period of Pm — 101 min. The amplitude and 
complexity of the stray light background variations depend on the season of observing, 
the location of the target star relative to the bright limb of the Earth and the orientation 
(roll) of the spacecraft. In the case of the V471 Tau photometry, the periodic fluctuations in 
background translated into photometric uncertainties in the stellar signal ranging from point- 
to-point scatter with cr ~ 0.003 (about 3 mmag) at stray light minimum to a point-to-point 
scatter of about cr > 0.1 at stray light maximum. 

The dark and flat field corrections were performed by monitoring individual pixel re- 
sponses during test exposures on fields empty of stars bright enough to rise above the back- 
ground. Photometry was extracted from the stellar images using a Moffat-profile point 
spread function model (IMoffat Ill969l ). The correlation in the raw photometry between the 
instrumen tal magnitude ligh t curve and the estimated sky background was removed as de- 
scribed in iRowe et al.l (120061 ). About 29% of the total number of data points were rejected 
because of pixel saturation during phases of the highest stray light in the MOST orbit and 
high cosmic ray fluxes when MOST passed through the South Atlantic Anomaly, as indicated 
by the orbital model of the local magnetic field strength. Additionally, about 6% of data 
points were rejected because of the relative uncertainty exceeding a = 0.015 of the mean 
light level. 
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The reduction and selection procedure left 56,383 measurements containing gaps of 
variable length spaced roughly by the MOST orbital period, resulting in a net duty cycle of 
about 65%. (We later conducted a period search after an even stricter selection of the data, 
with a duty cycle of 59%, as described in Section I3.3[ ) The time sampling and duty cycle 
provide excellent coverage in binary orbital phase and during the eclipses of V471 Tau. Note 
that the orbital period of the binary of close to 1/2 day always created a phase-coverage 
problem for ground based observations; the MOST data are entirely free of this limitation. 
The MOST photometry data (raw, and the reduced light curve used for analysis in this 
paper) are available in the MOST Public Data Archive on the Science page of the MOST 
web site: www.astro.ubc.ca/MOST. 



2.2. V471 Tau light curve 



The 19 orbital cycles of the binary monitored by MOST allowed us to investigate changes 
in the light curve from cycle to cycle, which is normally int erpreted as migr ation and evolution 
of spots on the magnetically active K dwarf component fjlbanoglulll978l ). The MOST data 
were divided into 1-day long subsets and those subsets were phased with the known orbital 
period of V471 Tau. Three of these subsets, from the beginning, middle and end of the 
10-day run, are overplotted in Figure [TJ A subtle trend is visible in that the rising portion 
of the light curve (in the phase interval 0.05 - 0.25) moves systematically later in phase with 
time, by a total of about 0.04 over 10 days. There is some evidence of this shift during the 
falling portion of the curve in the phase interval 0.75 - 0.95, but it is less pronounced. No 
phase shift is seen in the phase range 0.3 - 0.7, within the photometric scatter. 

The change s seen in the MO ST photometry resemble the "photometric wave migration" 
first reported by llbanoglul (119781 ) and discussed below in Section [3711 The average shift of the 
wave of ~ 0.002 phase/day indicates that it would take 500 ±250 days for the wave to make 
a full revolution {Pmigr)- This is somewhat longe r than the wave migration period f ound in 
previous studies (from ~ 180 d by llbanoglul (Il989l ) to 372 d by lSkillman &: PattersonI (Il988l )). 
although our estimate of the drift rate, based on only 19 orbital cycles, is necessarily crude. 

Seeing that the systematic changes in the light curve during the 10-day span of our 
observations were relatively small, with apparent shifts less than 0.01 mag at a given orbital 
phase, we calculated a mean light curve from the entire time series. This is presented in 
Figure [2l 



- 6 - 



2.3. DDO spectroscopy 

We obtained ground based spectroscopy of V471 Tau which partially overlapped with 
the MOST photometric run during 7-19 December 2005 (see Table [1]). A total of 37 
spectra at a dispersion of 0.14 A/pixel were collected using the Cassegrain spectrograph 
of the 1.88-m telescope at the David Dunlap Observatory. Since we expected the K-type 
dwarf in the system to dominate the flux at optical wavelengths, the wavelength range of 
the spectra was centered at Ha line, covering a red spectral window between 6425 and 
6705 A, (Figure [3]). This region contains a multitude of telluric lines which were removed 
during standard reductions performed using IRAF@ routines. The spectra were taken with 
the integration times and at intervals of about 30 minutes and could not cover all orbital 
phases of the binary because of the night/day breaks, commensurability of the binary period 
with one day and interruptions due to weather. The long integration times preclude any use 
of the spectroscopic data for improvement of the eclipse timing described in Section 13.21 



INTERPRETING THE LIGHT CURVE 



The phase diagram of the mea n light curve present ed in Figure [2] was modeled us- 
ing the PHOEBE software package (jPrsa fc Zwitterll2005l ). based on the Wilson- Devinney 
model. The orbital a nd physical parameters of both stars in the system were adopted from 
O'Brien et all J200lh : Rk = 0.96 i?©, Tk = 5,0A0K, Rwd = 0.0107 i?0. Two = 34,500 7^, 



a = 3.30 i?Q, i = 77.4°; the subscripts K and WD signify the K and white dwarf compo- 
nents, respectively. The atmospheric parameters for the red dwarf component were set to 
typical values for a K dwarf; limb darkening = 0.52, gravity darkening = 0.32 and albedo = 
0.5. 

The resulting model reproduces the general nature and amplitudes of the double-wave 
variability, and the depth of the eclipse, seen in the MOST light curve, as shown in Figure [2J 
It consists of the dominant smooth, wave-like variability and a relatively shallow (0.022-mag 
deep) total eclipse lasting 46.9 minutes, with steep shoulders each approximately 50 sec 
long. The photometric double wave is caused by ellipsoidal distortion of the K dwarf, with a 
minor modification due to the reflection effect. The asymmetry in the ellipsoidal distortion 
variability is believed to be due to spots on the K dwarf. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



- 7- 



3.1. Spot coverage 



In order to study the locations and extent of spots on the surface of the K dwarf, we 
used the residuals between the observed light curve and the modeled light curve (Figure H]) 
to estimate the required changes of the spot filling factor with orbital phase. Because of the 
orbital inclination of 77°, it is only possible to estimate changes in the mean spot coverage 
on the K dwar f disk within t h e lati tude range of —77° to +77°. Recent Doppler imaging 
observations of iHussain et al.l (120061 ) revealed that the K dwarf is rotating rigidly; this is 
confirmed by our determination of Vft-sini (Section 14.21) . As our run duration w as only 2.5 
times longer than the time span of the observation used by iHussain et al.l (120061 ) , we expect 
any changes of filling factor at a given phase to refiect spot rearrangement caused by the star 
activity rather than the star differential rotation. Also, any spot interpretation can address 
only the part of the spot coverage which is longitudinally asymmetric. 

Our results indicate that the smallest spot coverage occurred during the orbital phases 
0.6 — 0.7, while the largest occurred during 0.2 — 0.3. We seem to have observed a totally 
diff erent level of act i vity in the K dwarf than s e en du ring the Doppler imaging observations 
by iRamseyer et al.l (119951 ) and IHussain et al.l (120061 ). Our estimate of the peak-to-peak 
amplitude of the spot filling factor, 0.02 — 0.03 (depending on the assumed spot temperature 
differential values of AT = 2,000 — l, 000i^, as shown in F igure HI), is many times smaller 
than the changes of ~ 0.15 observed by iHussain et al.l (120061 ) in November 200 2. Also, in our 



data, the maximum spot coverage is inferred close to orbital phase 0.25, while IHussain et al. 



(120061 ) found the maximum around phase 0.07. The evolution of the spot coverage during 
the 10-day MOST observing run was still smaller, typically at a level of < 0.01, depending 
on the phase. 



3.2. Eclipse timing 



The relatively slow 10-sec photometric sampling rate (compared with the eclipse branch 
duration of 50 sec) and the temporal gaps left after the data selection made it impossible 
to measure times of individual eclipses accurate to a few seconds. Instead, we calculated 
the average eclipse time on the basis of a phased light curve of the entire time series to 
compare with earlier eclip se times in the li t eratur e. The phases were computed with the 
linear ephemeris given by iGuinan fc Ribad (120011 ). Because previous eclipses have been 
observed over a long time span (~ 35 yr) and the orbital period of the binary is short, we 
adopted a uniform time s ystem o f Helio centric Julian Date based on the Terrestrial Time 
(HJED), as advocated by iBastianI (120001 ). 
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The echpse time was determined after correction of the hght curve for the local slope 
created by the photometric wave. Since all contacts of the eclipse are not well defined (see 
Figure we determined the intersections of the averaged eclipse branches with horizontal 
line at the mid-depth level. The mid-p oint of both inter s ection s corresponds to the mid-point 
of the eclipse. With the ephemeris of iGuinan fc Ribad (120011 ). our mean epoch corresponds 
to = 25, 135. The shift in the observed time of t he mi d-point of eclipse is large compared 
to the predicted zero phase by iGuinan fc Ribad (120011 ): O — C = +248 ± 7 seconds, or 
over 4 minutes (see Figure E]). The MOST eclip se time determ i nation is shown compared 
to all available published data (as discussed by Ibanoglu et al. ( 20051 )) in Figure [61 The 
O — C curve continues an upward trend seen for about the last 10,000 orbital cycles. The 
implications of the MOST timing point are explored below. 



3.2.1. Third body 



The V471 Tau period changes visible in the eclipse — C diagram have been interpreted 
previously by several others as a light-travel-time effect caused by a perturbing third star in a 
long- period orbit in the system (jibanoglu et al.lll994j : iGuinan &: Ribaal200ll : llbanoglu et al. 



20051 ). This explanation is attractive because it could be reconciled with the main features 
of the — C diagram. It is also exciting because the mass of the hypothetical third body 
would be sub-stellar for a large range of possible inclination angles. 

Our new eclipse timing measurement shows that the long-anticipated downward bend 
in the — C diagram has not yet happ ened. Moreover, it d eviates substantially from the 
most recent third-body model proposed (llbanoglu et al.l 120051 ) by 52 sec, which is 3.6 times 
larger than a of the residuals for this model, as shown in the lower panel of Figure O Indeed, 
the MOST point is the largest deviation from this model so far observed. 

Therefore, we decided to recalculate the third-body model utilizing the same formalism 



as m llbanoglu et al.l (120051 ). With the new data augmented by the MOST result, the semi- 
amplitude of the — C variations, the third-body orbital period and its mass function are 
all slightly larger and longer than those given by llbanoglu et al.l (120051 ): see Table [2] for 
the full set of fitted parameters. The third-body orbital fit, although formally appropriate, 
remains uncertain because we still do not see the bend in the — C curve. In fact, as 
is shown in Subsection 13.2.31 below, it is reasonable to assume that the period has been 
constant since E ^ 15,000, i.e., over the last ~ 14 years. However, if we continue to see 
a straight line in future extensions of the — C diagram, this will not necessarily exclude 
the third-body hypothesis. Figure [6] includes a fit to a third-body model whose orbit has an 
even longer period which can still match the observations. Note that the orbital inclination 
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range necessary to preserve the sub-stellar mass of the third body will decrease to a very 
small range of angles if the current linear trend in the O — C deviations continues. 

The suggested parameters of the hypothetical third body in the V471 Tau system indi- 
cate that this object may be detectable with modern infrared teles copes or interferometers . 



With a larger mass function and a longer orbital period than in iGuinan fc Ribad (120011 ). 
the separation and brightness of the third body can be as large as 0.9 arc second and K ~ 
13.3 mag; see Table E] for predictions of the third body parameters for its different orbital 
inclinations. 



3.2.2. Apsidal motion 



If the binary orbit is even slight ly ellipso i dal, it may show a motion of the l ir ie of apses. 
This explanation was mentioned by iHerczegl (119751 ) and ISkillman fc Patterson! (119881 ). but 
then dismissed as an unlikely cause for the changes in eclipse times. We performed a least- 
squares fit of the O — C curve with the first-order formula given by iTodoranI (119721 ) and 
found that a very narrow range of eccentricity, e = 0.0121 ± 0.0006 (with 98% confidence 
level) , is required to explain the latest O — C results we have presented. See Table |3] and 
Figure [71 Although the orbit is expected to circularize in a very close binary system like 
V471 Tau, our fit to a slightly non-zero eccentricity is surprisingly close to the one we find 
from our radial velocity orbital measurements (see Section WA\ below) . 



3.2.3. Sudden period changes 

Without assuming anything about the actual nature of the O — C changes, it may be 
argued that the curve is composed of a few straight-line segments, each corresponding to a 
constant period, and of relatively short intervals where abrupt period changes take place. 
The portions of the O — C diagram from epochs E 2, 500 to 10, 500 and from E ^ 15, 000 
onwards appear to be consistent with two different constant periods. Least-squares linear 
fits to both segments of the O - C diagram yield periods of 0.52118305(4) and 0.52118384(4) 
days, respectively (the errors in parentheses are given in units of the last decimal place), 
corresponding to a relative period change of AP/P ~ 1.5 x 10~^. 

A sudden period change may be explained as a result of mass transfer or mass loss 
in a binary. For V471 Tau, we do not know if the possible donor, the K dwarf, is more 



massive than the mass recipient, the WD, but this is the most probable (jO'Brien et al. 



200ll ). In that case, the favorable scenario of a recent period increase is mass loss at the 
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level of ~ 3.8 x^OlLMa/w- flHilditchl 120011) . Taking the masses of both components at the 
limits of the lO'Brien et al.l (120011 ) ranges we can also consider the case when the donor is 
the less massive star. Such a situation would require conservative mass transfer at a level 
of ~ 3.6 X 10~^MQ/yr to explain the recent period increase. Both mass-loss rates appear 
to be large and unlikely for V471 Tau as they would result in other detectable phenomena. 
Moreover, both a period increase and a period decrease have been observed for the system 
so the complete picture would have to be even more complex. 

The latest period change took place over some AE ~ 2500 cycles so the inferred time 
scale, T, was T = {dlnP/dt)~^ ~ 2 x 10^ years. This is a relatively short time scale for any 
thermal equilibrium adjustment in the K dwarf, but of course may relate only to the outer 
layers of its envelope. 

The standard deviation in the residuals of the second segment of a = 22.7 s (Figure [7]) 
is slightly larger than for any of the previous fits (14.9 s for the third-body model and 
16.6 s for the apsidal motion model) but is still acceptable if superimposed upon possible 
short-timescale variations which are considered below. 



3.2.4- Periodic residuals from eclipse timing models 

Every one of the O — C models we calculated generates residuals with a larger than 
the accuracy of the eclipse timings (typically a few seconds). We performed a search for 
periodicities in the residuals and found that regardless of the model used, there is evidence 
for a 10-year period in the timing residuals. 

To investigate this further, we decided to employ a multi-harrnonic analysis of variance 



(MAOV) with 2 harmonics, as described in ISchwarzenberg-Czernyl (119961 ). This method uses 



orthogonal polynomials to model the data and the analysis of variance statistics to evaluate 
the quality of the result. The MAOV periodogram obeys Fisher's probability distribution 
with 2A^ + 1 and K — 2N — 1 degrees of freedom, where N is the number of harmonics used 
and K is the number of observations. The quantity F{2N+ 1, K — 2N — 1) shown in Figure[8] 
measures the ratio of powers of the periodic signal and residual, fully random noise. 

The amplitude of the variations we find in the O — C residuals is similar for all three 
models we adopted, at the level of 20-25 s and indeed indicates a ty pical underlying variation 



with a time scale of about 10 years. The 5.5-yr period found by llbanoglu et al.l (120051 ) - 
which was connected with the ~5-yr period in the mean brightness variations of the system 
- is also present, but at a much lower significance level (see Figure [8]). 
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3.3. Short-period oscillations 

Fluctuations with a period of 555 s were discovere d in soft X-ray flux from the V471 Tau 
system by the EXOSAT satellite (I Jensen et al.lll986l ). In 1991, 131 h ours of continuous U- 
band photometry of V471 Tau by the Whole Earth Telescope (WET) (jClemens et al.lll992l ) 
resulted in the detection of three periods: 554.63, 561.59 and 277.319 s. The dominant 555-s 
variability (with its 277s harmonic) was attributed directly to the accreting magnetic polar 
caps on the white dwarf component of the system, and the 562s signal to the same radiation 
reprocessed by the K dwarf atmosphere. 

To search for short-period variations in the MOST photometry, we first removed varia- 
tions caused by the binary revolution and rotation of the spotted component. The data were 
"rectified" by fitting the data with least-squares low-order polynomials and then dividing 
by the fitted function. The eclipses and flare events (see Section [3^ below) . accounting for 
about 7% of the total time series, were excluded from the fit, resulting in a net duty-cycle of 
59%. The remaining 52,371 brightness measurements of the binary, as well as correspond- 
ing measurements of both reference stars, were used to calculate MAOV periodograms, as 
described above in Subsection I3.2.4[ 

Analysis of the resulting periodogram revealed that none of the three WET periods is 
present in the MOST data, but their absence is easy to understand. While the white dwarf 
contribution to the total brightness of the system in the U band is about 39%, it is only 2.3% 
in the broad MOST photometric bandpass which includes considerable red flux. Therefore, 
the relative amplitude of the variations in MOST photometry is expected to be about 17 
times smaller than in WET photometry. The relative signal would be ~1.8 x 10~^, which 
is slightly below our estimated, one-sigma detection limit of about 2 x 10^^ in these data. 
This value was calculated by folding the data with a period incommensurate with any of the 
V471 Tau variations and MOST orbital harmonics. The noise estimation was also confirmed 
with the photometric data of both reference stars. 

Thus, the non-detection of the white dwarf pulsations in the broad MOST passband is 
entirely predictable. We can conclude only that the pulse amplitude (and presumably the 
polar accretion rate) did not increase significantly since the WET campaign in 1991. 



3.4. Flare activity 



Se veral flare-like events h ave been repor ted in V471 Tau by^ucinski ( 1981 ). Tunca et al. 



( 119931 ) ■ llbanoglu et al.l (120051 ) and others. lYoung et al. Jl983[ ) bund that flares are most 
likely to occur when the brightness of the system is near its minimum, when the K dwarf 
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was thought to have its most spotted hemisphere facing Earth. 

In the MOST hght curve, we identified seven events we would consider fiare-hke, al- 
though two of them were only partially recorded due to gaps in the data. This is the first 
detection of white-light flares by the MOST satellite and probably the largest homogeneous 
set of V471 Tau flare-like events observed so far. The durations of these events varied from 
about 10 to over 35 minutes, but their shapes all share the same rapid rise and slower decay 
characteristic of flares seen in visible light. The candidate events are shown in Figured 



In contrast to lYoung et al.l (Il983l ). we did not flnd any correlation of the flare events 
with the photometric wave minimum. The flares occurred during phases of the lowest as 
well as the highest spottedness of the K dwarf, with no apparent concentration in phase. 
The symbols at the bottom of Figure H] mark the phas es when the fl a res o ccurred. Using 



luminosities of both components in the V band given by lO'Brien et al.l (l200ll ). we estimated 
a lower limit to the energy released during the whole duration of a typical flare observed 
during the MOST run at about 10^^ erg (see Table H]). The energies of each of the seven flares 



we ob served are comparable to the energy released by the flare reported by llbanoglu et al. 



(120051 ) and are at the top of the range of energies released by all flare-like events reported for 
V471 Tau. Because the activity cycle of V471 Tau still remains to be characterized in terms 
of its period and intensity, we cannot relate the observed incidence of flares to the phase in 
this cycle. We note only that all the observed flares share the shape, duration and energy 
with those reported for typical RS CVn systems. 

The number of detected flare-like events corresponds to a total number of about 10 such 
events during the 10-day span of the MOST observations. Considering the limitations of 
ground-based observations one would expect to be able to detect a maximum of 4 flare-like 
events during the same period of time. 



4. INTERPRETING THE SPECTRA 

The typical S/N of the DDO spectra of V471 Tau is about 30. The contribution of 
the white dwarf component to the total light in the observed wavelength range is less than 
1%, so its contributions to the spectroscopic analyses described below are negligible. Our 
discussion of the spectroscopic results is limited to the K dwarf in the system. 
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4.1. Radial velocities 



To derive t he radial veloci ties (RV) of the K dwarf, we used the Broadening Function 
(BF) technique (lRucinskilll999l ). Spectra of four different K-type standard stars (HD 62509, 
HD 65583, HD 3765, HD 103095) were adopted as templates. The resulting broadening 
functions were fitted by a rotational line-broadening profile, with a linear limb-darkening 
coefficient of 0.5 2 (assumed to be t ypical for a K-type dwarf in the observed wavelength 
range), following Ivan Hammd (119931 ). The resulting RV measurements are listed in Table [1] 



We performed two independent least-squares fits to the radial velociti es, assuming first 
a circu lar and then an eccentric orbit, at a fixed orbital period as given by lGuinan fc Ribas 
(I2OOII ). but with the time of conjunction taken from the MOST light curve. The results of 
the fits and their residuals are plotted in Figure [TOl The quality of both fits, evaluated by 
calculating the standard deviations of the residuals, is essentially identical for both types 
of orbits, with a ~ 1.25 km s~^. The fact that a is not reduced for a model with more 



free p arameters suggests that the eccentric orbit solution is not necessary (ILucy fc Sweeney 
I971I ). although obviously this is not a proof for perfect circularity of the V471 Tau orbit. 



All our orbit a l mod el parameters (Table [5]) agree very well with those obtained recently 
by iHussain et al.l (120061 ) . but t hey deviate slight ly from those obtained previously with the 
same DDO 1.88-m telescope by lBois et al.l (119881 ). The amplitude we find is larger by about 
1.5 — 2 km s~^, and the center-of-mass radial velocity is about 2 km s~^ smaller. 



4.2. Projected rotation velocity 



A bonus of the BF analysis is the availability of the projected rotation profile of the star 
onto radial velocity space (Figure [TT]). This shape can be interpreted through a solid-body 
rotation to estimate the projected equatorial velocity V/i-sini. In the BF determination, 
we used HD 3765 as a standard star because its spectral type, K2V, is identical to that 
of the V471 Tau K dwarf. An average of the projected rotational velocities for all spectra 
is V/^fSinz = 91.9 ± 2.5 km s~^. The value is corrected for the broadening introduced by 
the method, the magnitude of which can be estimated by applying the B F method to the 



templ ate i tself. The result is c onsistent with previous estimates made by iRamseyer et al. 



(I1995I ) and lHussain et al.l (120061 ) (91 ±4 and 91 ±2 km s ^, respectively) and all are consistent 
with synchronous rotation of the K dwarf in V471 Tau. 
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4.3. Ha emission 



The Ha line was det ected in emission i n V471 Tau by iLanning &: Etzell ( 19761) . Sub- 
sequent detailed studies ( Young et al. 1988 : Bois et al. 199irRottler et al. 2002 ) revealed 
orbital phase-dependence of the emission as well as long-term changes of its equivalent width. 

We extracted the Ha emission from the absorption profiles of our spectra by again us- 
ing the standard star HD 3765 as a template. HD 3765 has the same spectral ty pe as the 
V471 Tau K dwarf and rotates very slowly at Vsini ~ 1 km s~^ (ISoderblomlll985l ) . We con- 
volved the standard spectrum with the rotational profile calculated for V/^ sini = 92 km s~^ 
(our value for V471 Tau) and fitted the resulting modified spectrum to each of our V471 Tau 
spectra in two wavelength ranges: 6540 — 6555 A and 6570 — 6585 A (see Figures [3] and [T2|) . 
Subsequently, we used the net Ha emission to derive the radial velocities and equivalent 
widths of the emission line (Table [1]). The extracted Ha profiles were symmetrical thus 
allowing us to use a Gaussian fit for measuring RV and numerical integration for equivalent 
widths. The radial velocity of the Ha emission (Figure [T^ follows the K dwarf orbital vari- 
ations, but with a smaller amplitude of about 120 km s ~^, as estimated from a sinusoidal 
fit. Such behavior was observed during 1975 - 1984 by iBois et al.l (Il99ll ). but with a still 
much smaller amplitude of ~ 75 km s^^. We observe that the Ha equivalent width changes 
symmetrically with respect to its maximum at orb ital phase 0.5 (Figure [T51) . in a very simi- 
lar way to what was reported by iBois et al.l (Il99ll ). The amplitude of the equivalent width 
variability in our data is about 1.2 A with the maximum emission of about —0.5 A at phase 
0.5. 



Long-term changes of Ha emission were detected by iBois et al.l (|l99ll ). who showed 
that the emission strength diminis hed between 1 9 75 an d 1983 and then grew rapidly in 
1984. More recent observations by iRottler et al.l (120021 ) have shown that since 1985, the 
emission was dropping again, until it finally vanished in 1992. This suggests that the long- 
term variation in Ha emission strength may be periodic, with a period of roughly 9 years. 
Our measurements show that in December 2005, the emission strength was comparable to 
its average values in the past. This is consistent with a 9- year periodicity, since our DDO 
spectra were obtained about 2 years after the latest expected emission maximum in such a 
cycle. 



5. SUMMARY 



The nearly continuous MOST spacebased photometry of V471 Tau covering 10 days in 
December 2005, combined with partly simultaneous DDO groundbased spectroscopy, moni- 
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tored a fairly quiescent stage in the activity of the K dwarf in this close binary system. This 
is apparent in the light curve which deviates relatively little from the model and almost does 
not change during the whole observing run. Even during such a stable time, seven candidate 
flare events were observed in 10 days, whose estimated energies would be among the highest 
ever seen in V471 Tau. There is no correlation between the times of the flares and orbital 
phase. 

The main features of the orbital phase diagram of the MOST photometry are well 
reproduced by our eclipsing binary light curve synthesis model. The largest systematic 
deviation in the double-wave light curve is only about 0.02 — 0.03 mag and is consistent 
with spots on the K dwarf which is expected to rotate synchronously with the orbit. The 
amount of spottedness on the star did not change much during the MOST observing run - 
by no more than about 1%. This supports the claim that the K dwarf was observed close 
to a minimum in its activity cyc le. A half-orbi t al-per iod modulation of the radial velocity 



residuals was reported earlier by iHussain et al.l (120061 ) and interpreted as an asymmetry in 



spot distribution on the K star's surface. We see no such resid u als in our radial velocity 



measurements. We note that the residuals seen by lHussain et al.l (l2006l ). the radial velocity 
curve we obtain, and the O — C variations in echpse times observed over the past 35 years, 
could all be interpreted as a small non-zero eccentricity of the orbit of V471 Tau. 

Because of the broad bandpass of the MOST photometry with substantial flux in the 
red, and the red wavelength range of the DDO spectra, the white dwarf contributes only 
about 2% and 1% of the total intensity of the system, respectively. We were therefore unable 
to constrain the properties of t he hot white dwarf in the system or confirm the oscillation 



frequencies detected by WET ( IClemens et al.l Il992l ) , since the relative amplitudes in the 



custom-filter, broadband MOST photometry would be about 17 times smaller than in the 
WET tZ-band photometry. The positive aspect of this is that our estimates of the K dwarf 
properties from MOST photometry and DDO spectroscopy are not contaminated by the 
white dwarf, but we can use the timing of the white dwarf eclipses to investigate aspects of 
the orbit of the V471 Tau system. 

Changes in the O — C values of the times of eclipse of the white dwarf can, however, be 
explained by at least three entirely different models: (1) There could have been at least two 
abrupt period changes in the orbit of the system in the last 35 years, although there is no 
obvious mechanism for this. (2) There could be apsidal motion due to a slightly eccentric 
orbit. (3) The V471 Tau system might be a trinary, with a third low-mass companion in 
a long-period orbit. The last two periodic phenomena both predict that the O — C eclipse 
timing deviations must drop in the future (see Figure |7]). The small eccentricity which could 
explain the — C diagram is also in agreement with the formal solution of the radial velocity 
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curve of the K dwarf from our high-quality DDO spectra, but its value is currently below the 
direct spectroscopic detection threshold. Future accurate eclipse timing observations, such 
as performed by the MOST satellite, are desired as they may resolve the dilemma between 
those three models. 

The O — C residuals do show a convincing residual periodic variation with a period 
of about 10 years, regardless of the model used to explain the longer-term changes. This 
variation may be due to an activity cycle in the K dwarf, but this is a highly speculative 
explanation. We note that the Ha emission appears to change in intensity in a characteristic 
time scale of about 9 years, perhaps coincident with the periodicity in eclipse time variations 
at the frequency resolution of the entire data sample at hand. The 10-year period in the 
O — C residuals ma y also be related with the 5. 5- year period in the system mean brightness 
variations found by llbanoglu et al.l (120051 ) as its multiple. Nevertheless we think that both 
periods are too uncertain to firmly connect them at this stage of the study of V471 Tau. 

The DDO spectra yield a new radial velocity curve for the orbit of the K dwarf, and an 
improved determination of the projected rotation of the star, V^sini = 92 km s~^ based on 
high-quality BF (broadening function) profiles. The spectra also enabled us to measure the 
Ha emission velocities and changes in its equivalent width. The H a emission of V471 Tau 
showed the s a me o rbital phase dependence as observed before by iBois et al.l (Il99ll ) and 
Rottler et al.l (|2002| ) with maximum emission at phase ~0.5. The observed amplitude of 
equivalent width variations of about 1.2 A was average for the system and consistent with 
the 9-year period noted by previous investigators. Unfortunately, the 13- year gap between the 
most recent published Ha emission observations of V471 Tau and our new DDO observations 
does not allow us to reliably verify the periodic character of the mean emission strength 
variation. 

A new feature of the Ha emission revealed by our observations was its much larger 
amplitude of radial velo city variati o n (12 km s~^) compared to that reported by earlier 
observers (75 km s~^ by iBois et al.l (Il99ll )). This suggests that the source of the emission 
was less concentrated around the sub-white-dwarf point on the K star as had been seen in 
the previous data. 
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Table 1. Spectroscopic observations. 



No. 


H JED -2,453,700 


photometric 


S/N 




Ha EW 


Vrad of Ha 






phase 




(km s-l) 


(A) 


net emission (km s~^) 


1 


17.00275 


0.5162 


30 


16.8 


-0.578 


42 


2 


17.715D1 


0.7717 


30 


-113.2 


0.183 


-77 


3 


17.72do2 


{J.7922 


30 


-108.3 


0.249 


-57 


4 


17.7-3841 


0.8154 


30 


-101.1 


0.361 


-77 


5 


1 / . /4yoo 


0.8364 


30 


-91.2 


0.407 


-70 


6 


17.7d19d 


0.8606 


30 


-78.8 


0.493 


-22 


7 


17.77283 


0.8815 


30 


-66.3 


0.620 


-7 


8 


1 T TO /loo 

17.7o4o2 


0.9045 


12 


-50.3 


0.801 


-11 


9 


17.79054 


0.9250 


30 


-31.3 


0.727 


-7 


10 


17.8U788 


0.9487 


30 


-12.4 


0.789 


36 


11 


17.81871 


0.9695 


8 


12.2 


0.745 


15 


12 


17.83145 


0.9939 


8 


31.0 


0.538 


57 


13 




0.2443 


20 


184.8 


0.165 


167 


14 


18.49374 


0.2647 


30 


184.2 


0.103 


186 


15 


18.5U572 


0.2877 


30 


183.2 


0.025 


168 


16 


18.51648 


0.3083 


30 


176.1 


-0.050 


166 


17 


18.52928 


0.3329 


30 


165.2 


-0.143 


159 


18 


18.54004 


0.3535 


30 


157.4 


-0.207 


141 


19 


18.55221 


0.3768 


30 


142.2 


-0.278 


147 


20 


18.56299 


0.3975 


30 


127.6 


-0.375 


129 


21 


18.57489 


0.4203 


30 


107.5 


-0.411 


110 


22 


18.58558 


0.4409 


30 


91.8 


-0.485 


92 


23 


18.59955 


0.4677 


30 


63.3 


-0.535 


81 


24 


18.61024 


0.4882 


30 


44.4 


-0.544 


62 


25 


18.62229 


0.5113 


30 


22.4 


-0.485 


42 


26 


18.63295 


0.5318 


30 


4.8 


-0.520 


26 


27 


18.64494 


0.5548 


30 


-16.4 


-0.388 


18 


28 


18.65570 


0.5754 


30 


-34.7 


-0.319 


3 


29 


18.66658 


0.5976 


15 


-52.8 


-0.298 


-22 


30 


18.67986 


0.6218 


30 


-70.3 


-0.135 


-25 


31 


18.69057 


0.6423 


30 


-84.7 


-0.102 


-63 


32 


18.70082 


0.6635 


30 


-92.5 


-0.121 


-42 


33 


23.79578 


0.4377 


30 


91.5 


-0.393 


100 


34 


23.80649 


0.4583 


30 


75.0 


-0.447 


90 


35 


23.81862 


0.4816 


10 


53.2 


-0.394 


61 


36 


23.82925 


0.5020 


10 


32.9 


-0.440 


55 


37 


23.84138 


0.5252 


20 


6.5 


-0.403 


35 
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Table 2. The best fit parameters for a tliird body model. 



parameter 


value 


To (HJED) 


2440610.06446 ± 0.00008 


^'O (d) 


0.521183449 ± 0.000000008 


Ps (yr) 


33.7 ± 0.9 


0.12 sin(i3) (AU) 


0.32 ± 0.02 


semiamplitude (sec) 


159 ± 6 


e 


0.28 ± 0.09 


u}3 (degr) 


79 ± 10 


/(Ms) (Mq) 


(2.9 ± 0.3) ■ 10-5 



Note. — The errors were estimated using the boot- 
strap sampling method for the 98% confidence level; 
this is why they appear t o be large when comp ared to 
other fits, for example bv llbanoglu et al] ll2005h . 



Table 3. The best fit parameters for the orbital apsidal motion. 



parameter value 



P (day) 


38.884 ± 


0.007 


semiamplitude (sec) 


173 ± 9 




e 


0.0121 ± 


0.0006 


UJO (degr) 


174 ± 3 




uJl (degr/day) 


0.025348 


± 0.000005 



Table 4. A list of flare-like events observed by MOST. 



nr 


start time (HJED) 


phase 


duration (min) 




Emin (erg) 


1 


2453709.164 


0.379 


19 


0.014 


9.8 • 10^3 


2 


2453709.351 


0.738 


12 


0.014 


3.2 • 1033 


3* 


2453710.59 


0.12 


> 16 


> 0.016 


> 1.8 ■ 1034 


4 


2453711.345 


0.564 


13 


0.011 


1.0 ■ 1034 


5 


2453714.770 


0.135 


36 


0.018 


2.9 • 1034 


6* 


2453717.01 


0.43 


> 33 


0.018 


> 4.4 ■ 1034 


7 


2453718.156 


0.632 


10 


0.017 


1.9 • 1034 



Note. — Emin designates the lower limit for total energy released in the V band. 
The fiares marked with a star were observed only partially. 
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Table 5. Parameters of spectroscopic orbits. 



parameter 


circular orbit 


eccentric orbit 


Kk (km s-l) 


150.5 ±0.4 


150.2 ±0.5 


asini {Rq) 


1.550 ±0.004 


1.547 ±0.005 


e 




0.012 ±0.003 


u)\ (degr) 




75± 18 


Vo (km s-i) 


35.7 ±0.3 


35.2 ±0.3 


a (km s^-"^) 


1.25 


1.25 



Table 6. Parameters of hypothetical third body in the V471 Tau system. 



«3 (degr) 


Mz (Mq) 


T.ff (K) 


IoqL/Lq 


my 


ruK 


dmax (mas) 


Tmax (year) 


85 


0.045 


1540 


-4.3 


29.6 


15.4 


420 


2014.1 


60 


0.052 


1730 


-4.1 


26.8 


15.0 


490 


2015.0 


45 


0.064 


2060 


-3.8 


24.0 


14.4 


610 


2016.1 


30 


0.090 


2660 


-3.2 


19.2 


13.3 


910 


2019.0 



Note. — The ph ysical parameters of the third body are based on nongray dusty models of 
IChabrier et all 1 I2OOOI) . assuming the age of the system of 625 Myr and the dystance of 46.8 pc. dmax 
designates maximum apparent separation between the V471 Tau binary and the third component. 
Tmax is the time of the nearest maximum separation. 
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Fig. 1. — The mean V471 Tau light curve, averaged in phase with 0.01 phase bins, for 3 
selected days of the MOST observations at the beginning, middle and end of the run. A 
lack of obvious changes in the light curve during our observations beyond the global shifting 
at phases ~ 0.75 — 1.25, can be interpreted as a relatively low activity in the spot re- 
arrangement. 
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Fig. 2. — The mean V471 Tau light curve averaged in phase with 0. 001 phase bins (dots) . 
The hne shows the hght curve calculated by the PHOEBE package (jPrsa fc Zwitterl |2005| ) . 
based on the published V471 Tau parameters. Deviations caused by the spots on the K dwarf 
are not included in the model. 
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Fig. 3. — Comparison of the DDO spectra of V471 Tau. Top: The average of the 27 best 
quahty spectra after correction shifts for the orbital motion. The Ha hne was omitted due 
to its variabihty (see Subsection 14.31) Middle: The standard star HD 3765 spectrum after 
convolution with the broadening profile. Bottom: The HD 3765 spectrum as observed. 
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Fig. 4. — Relative changes of the spot-filhng factor versus the orbital phase for different spot 
temperatures, as indicated on the right vertical axis. The curve was obtained by comparing 
observed light curve (averaged in phase with 0.001 phase bins) with the theoretical one 
calculated with the PHOEBE package, as described in the text. The crosses on the bottom 
axis mark the phases of seven detected flare-like events (see Subsection 13. 4p . 
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Fig. 5. — All observations used for the V471 Tau eclipse timing (small dots) are shown in the 
phase diagram together with t he running averag e data (large dots). The phase shift of the 
mid-eclipse time relative to the iGuinan fc Ribad (120011 ) linear ephemeris is clearly visible. 
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Fig. 6. — The third-body model fits to eclipse timing observations of V471 Tau from the 
literature (filled circles). The open circle is the new timing from the MOST observations. 
A hypothetical, illustrative trend was created by adding a point of O — C = 0.004480 at 
E = 28, 000, that would follow the curve growing trend. The bottom plot show's the residuals 
of all available data with respect to the model predictions of llbanoglu et al.l ([20051). 
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Fig. 7. — Comparison of the third body, apsidal motion and straight line model fits to the 
available eclipse time observations of V471 Tau. 
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Fig. 8. — Multiharmonic analysis of variance periodograms (jSchwarzenberg-Czernyi Il996h . 
with frequencies up to 2 cycles per year, for different model residuals, as discussed in Sub- 
section 13. 2[ Horizontal lines show the levels of 0.001 (dashed line) and 0.01 (dotted line) 
probability of false detection. The most significant peaks appear around the same frequency 
~ 0.1 c/y for all models (the top two peri odograms also show it s ahas at ~ 0.05 c/y). Note 
that the 5.5 yr period (0.18 c/y) found by llbanoglu et al.l (120051 ) appears to be also present. 
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Fig. 9. — The seven flare-like events on V471 Tau which were detected during MOST obser- 
vation period. 
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Fig. 10. — The radial velocity curve for the K dwarf component of V471 Tau binary. The 
bottom panels show residuals for the circular and elliptical models, respectively. 
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Fig. 11. — The average broadening function of V471 Tau spectra derived with the standard 
velocity star HD 3765 of the same spectral type (sohd line). This BF profile was fitted by 
the rotational broadening profile to estimate the projected rotation velocity of the K dwarf 
component (dashed line). 
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Fig. 12. — A collection of representative spectra of V471 Tau taken at different phases. The 
variable strength and shifts in position of the Ha emission are clearly visible. 
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Fig. 13. — Variations of the Ha emission line. Top: The radial velocity changes of the net 
emission with phase, compared with those of the K-dwarf itself (the absorption spectrum). 
Note the reduced amplitude of about 120 km s~^. Bottom: Changes of the Ha hne equivalent 
width. The strongest emission is visible when the K-dwarf component is seen in the upper 
conjunction (when the face illuminated by the WD is directed to the observer around the 
orbital phase of 0.5). The emission is practically undetectable during the opposite phases. 



